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A modified version of the Plasma Beat-Wave Accelerator scheme is introduced and analyzed, 
which is based on autoresonant phase-locking of the nonlinear Langmuir wave to the slowly chirped 
beat frequency of the driving lasers via adiabatic passage through resonance. This new scheme is 
designed to overcome some of the well-known limitations of previous approaches, namely relativis- 
tic detuning and nonlinear modulation or other non-uniformity or non-stationarity in the driven 
Langmuir wave amplitude, and sensitivity to frequency mismatch due to measurement uncertainties 
and density fluctuations and inhomogeneities. As in previous schemes, modulational instabilities 
of the ionic background ultimately limit the useful interaction time, but nevertheless peak electric 
fields at or approaching the wave-breaking limit seem readily attainable. Compared to traditional 
approaches, the autoresonant scheme achieves larger accelerating electric fields for given laser inten- 
sity, or comparable fields for less laser power; the plasma wave excitation is much more robust to 
variations or uncertainties in plasma density; it is largely insensitive to the precise choice of chirp 
rate, provided only that chirping is sufficiently slow; and the quality and uniformity of the resulting 
plasma wave and its suitability for accelerator applications may be superior. In underdense plasmas, 
the total frequency shift required is only of the order of a few percent of the laser carrier frequency, 
and for possible experimental proofs-of-principle, the scheme might be implemented with relatively 
little additional modification to existing systems based on either solid-state amplifiers and Chirped 
Pulse Amplification techniques, or, with somewhat greater technological effort, using a CO2 or other 
gas laser system. 
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I. INTRODUCTION AND OVERVIEW 



The Plasma Beat- Wave Accelerator (PBWA) was first proposed by Tajima and Dawson (TD) |lj as 
an alternative to the short-pulse Laser Wake-Field Accelerator (LWFA), based on earlier analysis of 
beat- wave excitation as a plasma probe or as a mechanism for plasma heating 0, 0, 0] ■ Subsequently 
the PBWA concept has been studied extensively theoretically, numerically, and experimentally @, ||| 
SEIlllllllElliElHHEH- (Fora review, see In the original scheme, two laser 

fields co-propagating in an underdense plasma are detuned from each other by a frequency shift close to 
the electron plasma frequency, so that the modulated envelope resulting from the beating between the 
two laser fields can act ponderomotively on the plasma electrons to resonantly excite a large-amplitude, 
high-phase-velocity plasma wave suitable for particle acceleration. 

For a fixed beat frequency, performance of the PBWA is constrained by what is now known as the 
Rosenbluth-Liu (RL) limit, after the pioneering study in As the plasma wave amplitude grows, 
relativistic detuning effects eventually prevent further growth of the peak longitudinal electric field E z 
beyond a maximum value i£ RL , which, for any realistic laser and plasma parameters, lies below the cold, 
non-relativistic, one-dimensional (ID), wave-breaking limit Eq |2(ll21|: 

|^|<^^o(f^^) 1/3 <^o^= (1) 
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where c is the speed of light, m is the mass of the electron, e is the magnitude of its electric charge; ui p 
is the cold, linear electron plasma frequency, defined in Gaussian units as 

1/2 

uj p = (47m e 2 /TO) , (2) 

where n is the ambient electron number density; E\ and E 2 are the peak electric fields of the beating 
drive lasers; and ui\ and u>2 are their respective carrier frequencies. For the plasma waves of interest 
here, with high but sub-luminal phase velocities v p < c, the cold non-relativistic limit itself is smaller 
than the cold, relativistic (or nonlinear) wave- breaking limit [l^.l20j : 

E Q <E WB = [2{ lp -l)] 1/2 E , (3) 

which sets an upper bound on the amplitude of coherent plasma oscillations useful for particle acceler- 
ation in a cold plasma. Here 

7p =[l-^/c 2 ]- 1/2 (4) 

is the relativistic kinematic factor associated with the phase-velocity v p of the excited plasma wave, 
approximately equal to the characteristic group velocity v g (to be precisely defined below) associated 
with the envelope of the beating lasers, both of which will be nearly equal to the speed of light c 
for sufficiently underdense plasmas. (Large thermal spreads increase particle trapping and lower this 



tnermalsp: 

threshold for wave-breaking, as discussed in [^l^l^ .l 

It is not difficult to see the origin of this detuning effect. The effective nonlinear plasma frequency, in 
the presence of the laser drive and a plasma wave, is approximately given by 



'f + Sn/n ) 



7 PNL ~ p- W P) 

V /RMS 



(5) 



where Sn is an average density perturbation associated with the plasma wave, and 7 RM s is the root-mean- 
square relativistic factor of the electrons, including effects of both the transverse quiver in laser fields and 
the longitudinal velocity imparted by the Langmuir wave itself. Except for highly nonlinear situations 
where ponderomotive blow-out becomes important, the average density perturbation is typically small, 
and will be neglected here. The RMS transverse quiver velocity is roughly constant for the fixed- 
amplitude driving lasers typical of PBWA, but as the longitudinal motion of the perturbed electrons 
associated with the excited Langmuir wave becomes even weakly relativistic, 7 RM s will begin to increase, 
causing w PNL to decrease. Because of the sensitivity of resonance phenomena, 7 RM s need not grow much 
before this resonance shift severely limits the efficacy of the driving beat-wave. As this dynamical 
process is effectively reversible, if the detuning occurs before wave-breaking or other phase-mixing or 
dissipative effects, then the ponderomotively excited Langmuir wave is not actually saturated per se at 
some fixed value given by the RL limit, but exhibits a slow (compared to the plasma frequency) nonlinear 
modulation or beating of the wave amplitude, periodically peaking near the maximum value predicted 
by RL and then decreasing, as energy is exchanged between the Langmuir wave and the laser fields. 

Tang, Sprangle, and Sudan (TSS) |j],|8|] pointed out that one can, in principle, achieve somewhat higher 
amplitudes than that predicted by the RL limit by simply detuning the beat frequency to a value 
somewhere below the linear plasma frequency, matching the drive frequency to the nonline ar p lasma 
frequency for some non-trivial plasma wave amplitude. However, McKinstrie and Forslund [25j noted 
fundamental practical difficulties with this approach. For downward detunings at or beyond some small 
critical frequency shift, multiple solutions to the envelope equations appear, and the high-amplitude 
branch cannot be accessed reliably from the initially quiescent plasma state, especially in the presence of 
any small amount of damping. For sufficiently small detunings, moderate improvement over the RL limit 
is possible, but typically with peak amplitudes still lying below even the cold linear wave-breaking limit. 
Furthermore, as the detuning grows in magnitude toward this critical value, the excitation becomes 
increasingly non-monotonic, and the time needed to reach the peak amplitude increases and becomes 
comparable to the time-scales for ion modulational or decay instabilities that degrade the coherence of 
the plasma wave together with its utility for accelerator applications. 

Matte et al. (MMEBP) [HI suggested the use of a plasma with a time-varying density, carrying out the 
beat-wave excitation during the actual ionization process, when the free electron density is growing. The 
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increase in plasma density can then compensate for the increase in 7 RM s, and thus the effective plasma 
resonance may be maintained with a fixed driving beat frequency for a longer time. However, in order 
to yield appreciable improvement, such an approach would require impractically precise timing of lasers 
and some means to control the ionization rate, matching it, at least roughly, to the expected growth rate 
of the plasma wave. The usual method involving laser-induced field ionization leads to cross sections 
which are exponentially sensitive to the laser intensity and is therefore unsuitable for such an approach. 

Better still, Deutsch, Meerson, and Golub (DMG) [2l| suggested incorporating compensatory time- 
dependence into the drive lasers rather than in the plasma density. They proposed to partially overcome 
the relativistic detuning effect by using a chirped beat-wave excitation scheme, where the frequencies 
of one or more of the laser pulses are chirped downward starting from the linear resonance so as to 
compensate for the change in the nonlinear frequency of the growing plasma wave, ideally allowing 
amplitudes approaching the nonlinear wave-breaking limit. The nonlinear modulations are thereby 
reduced but not eliminated, consisting of a ringing about some non-zero saturation value rather than a 
full beating as in the original (RL/TD) scheme. 

Chirping too fast will cause the plasma oscillations to detune before they have a chance to grow much, 
while chirping too slowly appears to result in a slow, highly non-monotonic (ringing) excitation, where 
the long interaction times required to build up sufficient amplitude (on the order of a few ion plasma 
periods) can allow instabilities of the ionic background to grow and disrupt the coherence of the plasma 
wave US • Intuitively, one expects for the DMG approach that there exists an optimal intermediate 
chirp rate, or in fact even an entire optimal chirp profile, in which the driving beat frequency starts at 
the linear plasma resonance, where the initial coupling is strongest, and then self-consistently matches 
the instantaneous frequency to the detuning caused by the expected growth of the plasma wave. 

However, solving such a complicated nonlinear control problem so as to exactly match the nonlinear 
plasma frequency with the laser chirp, with the help of either numerical simulation or experimental 
calibration, would be both computationally daunting and experimentally impractical. The linear plasma 
frequency will be imperfectly known in practice, and even if it were obtained, the instantaneous value of 
the nonlinear plasma frequency actually depends on the initial absolute phase of the plasma wave, which 
is determined by the exact initial conditions of the plasma and the initial beat phase of the driving lasers, 
which are not normally controlled. However, DMG argue that such extrinsic frequency matching is not 
actually necessary, and that entrainment can be achieved without recourse to carefully-tailored pulses 
or to external feedback. They invoke the nonlinear dynamical phase-locking phenomenon known as 
autoresonance [2^. l30| , arguing that as long as the chirp is sufficiently slow (compared to the time-scales 
for the nonlinear growth and modulation described by RL) and in the right direction (namely, downward 
in frequency), then the externally forced, nonlinear dynamical system can self-adjust and automatically 
maintain approximately the desired phase entrainment between the driven wave and the ponderomotive 
drive, resulting in in an average increase in the oscillation amplitude. 

However, proposed as it was before the understanding of autoresonance had matured, especially with 
regard to the threshold behavior for establishing and maintaining entrainment and the importance of the 
initial detuning, the DMG scheme sometimes fails to produce appreciable phase-locking. This is largely 
due to its prescription of beginning the chirp on resonance, and becomes particularly problematic if the 
plasma density is imperfectly known or subject to shot-to-shot jitter, when it can lead to nonuniform 
growth or early saturation. In fact, there appears to be some uncertainty in the original description of 
the scheme |27| as to exactly how autoresonant the proposed excitation mechanism is, as to whether or 
to what extent "the precise form of the time-dependence of the frequency is inessential" or else there 
"exists an optimal chirp rate that provides the highest excitation rate." 

In fact the DMG chirped scheme can often achieve higher longitudinal electric fields than the original 
(RL/TD) approach, but does not do so universally or robustly, and cannot always produce fields ap- 
proaching even the linear wave-breaking limit for experimentally accessible parameter values. The peak 
amplitude and spatio-temporal uniformity of the final plasma wave may be sensitive to the exact chirp 
rate, to the initial detuning, and most especially to the precise value of the plasma frequency, which, 
once again, often is imperfectly known or subject to significant variation, due to limited measurement 
precision, to shot-to-shot jitter, or to single-shot density fluctuations. 

Informed by more recent results, we propose a novel variant of the chirped plasma PBWA concept 
which also exploits autoresonance, but enhanced via Adiabatic Passage Through Resonance (APTR) 
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|3ll l32t l3^|. Rather than chirping downward from the linear resonance, we actually start, counter to 
naive intuition, with a frequency shift well above the resonance, and then slowly sweep the beat frequency 
through and below resonance. With this approach, when starting from a quiescent plasma, the final 
state of the plasma wave is insensitive to the exact chirp history, and the excitation is more robust 
with respect to imprecise characterization of the plasma density, or actual shot-to-shot fluctuation of the 
density or single-shot variation in it, provided only that the spatio-temporal scales of variability are long 
compared to the frequency and wavenumber of the plasma wave, and the relative range of fluctuation 
or uncertainty is not too large. The chirp rate need not be matched a priori to any anticipated rate 
of growth in the plasma wave or resulting rate of relativistic detuning, but must only be chosen to 
be sufficiently slow so as to satisfy a certain adiabatic trapping condition and thereby ensure phase- 
locking. Because of the nonlinear nature of the ponderomotive forcing in the relativistic regime, this 
adiabaticity condition becomes increasingly difficult to maintain as the wave amplitude grows, and for 
a constant chirp rate eventually the plasma wave would fall out of phase-locking with the drive. This 
leads to a true saturation in amplitude, with significantly less of the undesirable ringing seen in the 
original RL/TD scheme. Before the onset of saturation, the frequency of the excited plasma wave can 
be closely entrained to the instantaneous beat frequency of the drivers, and the overall phase can remain 
reasonably entrained as well, while the plasma wave amplitude monotonically grows to automatically 
and self-consistently adjust itself to the monotonically decreasing beat frequency. 

Not only is our excitation scheme more robust, so too is our method of analysis. The Lagrangian fluid 
formalism developed by RL and then used by DMG employs a power series expansion which is valid 
only for weakly relativistic motion in both the transverse and longitudinal directions, thereby limiting 
its domain of applicability to laser fields of sufficiently low intensity (implying non-relativistic quiver 
motion), and to excited plasma waves of sufficiently low amplitude (with electric fields well below wave- 
breaking) . Therefore their equations of motion become increasingly untrustworthy precisely in the regime 
of interest, where the plasma wave amplitude grows to some appreciable fraction of the wave- breaking 
value Eq. Here, we instead employ a fully nonlinear, Eulerian fluid model that allows for arbitrarily 
relativistic electron motion below the nonlinear wave-breaking limit E WB , similar to that discussed in 
for general laser-plasma interactions and in f° r PBWA investigations, and also to a pioneering 
treatment in J^j. In the weakly relativistic limit, this Eulerian approach can be shown to agree exactly 
with the Lagrangian treatment of RL. 

This analytical fluid model, and the various physical assumptions which go into it, are discussed in 
Section [H] In order to treat the autoresonant nature of the problem more easily and transparently, we 
re-formulate the dynamical equations in a fully Hamiltonian form in Section IlIII In Section llVI we use 
this canonical formalism to analyze the autoresonant aspects of the beat-wave excitation, from the initial 
linear phase-locking regime through the weakly and strongly nonlinear trapped phases to non-adiabatic 
saturation. In Section[V] we discuss features of some realistic examples relevant for possible experimental 
implementation and investigation. Section IVII summarizes our preliminary assessment of the merits and 
limits of our autoresonant PBWA, especially in comparison to other beat-wave approaches. We then 
offer brief conclusions from our initial investigation and prospects for future study in Section IVIII 



II. FUNDAMENTAL EQUATIONS 

Our study of wake excitation is based on an approximate, but convenient and widely-used model. The 
plasma is treated as a cold, collisionless, fully relativistic electron fluid moving in a stationary, neutral- 
izing, ionic background, coupled to electromagnetic fields governed by Maxwell's equations. We restrict 
our analysis to one-dimensional (ID) geometry, where all dynamical quantities depend only on the lon- 
gitudinal position z and the time t, and we assume a completely homogeneous and initially quiescent 
background state of the plasma. The plasma is assumed to be highly underdense, i.e., ui p <C ui±, UJ2- 
In addition, we assume prescribed laser fields, neglecting throughout the entire interaction any changes 
to the laser envelopes due to linear effects such as diffraction or group-velocity dispersion, or any non- 
linear back-action of the plasma on the lasers such as depletion, self-focusing, photon acceleration due 
to ionization or density variation, as well as Raman scattering, self-modulation, and other instabilities 
[3{J|3(j. This model, although simplified, nevertheless reveals the essential features of autoresonance 
and its potential advantages for the PBWA. Possible extensions to more realistic models, as well as some 
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arguments for the validity of our general results beyond the strict applicability of this model, especially 
in the light of density fluctuations, are discussed in Sec. I VII and Sec. IVll| 

The cold, collisionless fluid treatment assumes that the electron temperature is sufficiently small so 
that: the thermal energy k^T e is negligible compared to the typical kinetic energy associated with the 
transverse quiver motion in the driving laser fields; thermal corrections to the Langmuir dispersion 
relation are small, which, for waves with relativistic phase velocities simply requires that k^,T e is much 
smaller than the electron rest energy mc 2 ; and collisional damping of the laser fields is small over the 
interaction time, as are Landau and collisional damping of the excited plasma wave. Neglect of ion 
dynamics is strictly valid provided the time-scale for ion motion, typically of the order of a few times 
wr , where 

Ui = (4Tm Z?e 2 /Mi) 1/2 (6) 

is the ion plasma frequency for ions of mass Mi and charge +Zie, remains longer than the duration T 
of the lasers: i.e., ujiT < 2tt. 

Within this model, the continuity, momentum, and Poisson equations, respectively, are, in Gaussian 
units: 

^l + V-(n e v) = (7) 

?T + (wV)p = efv*- - x (V x A)) (8) 

ot \ c / 

V 2 $ = -4we(n - n e ), (9) 

where A is the vector potential and $ the scalar potential in the Coulomb gauge; n e is the electron 
number density and no the background ion number density, assumed to be equal in the absence of 
perturbations; v is the electron velocity; and p = rwyv is the kinetic momentum, where 

7 = [i - M 2 /c 2 r i/2 (io) 

is the relativistic factor associated with electron motion. 

Neglect of diffraction or other transverse effects in the laser fields requires that characteristic laser spot 
size Wo sufficiently exceeds both the characteristic laser wavelength A and the transverse length-scale 
Ar± for variations in the plasma, while Ar± itself should greatly exceed the Langmuir wavelength 
A p ~ c/cjp to ensure the validity of our assumption of a ID, homogeneous medium. Neglect of the 
nonlinear evolution of the lasers effectively imposes certain constraints on the duration and intensity 
of the pulses, as summarized in 0. In the Coulomb gauge, the vector potential A is solenoidal (i.e., 
divergenceless) , which in ID also implies that it is geometrically transverse, so A = A±(z,t) will be 
polarized perpendicular to the propagation direction +z. For simplicity each drive laser will be assumed 
to possess right-handed circular polarization, and to have a flat-topped profile of total duration T, where 
the leading edge of each laser enters the left edge of the plasma at z = at the initial time t = Q. 

Defining a scaled, or dimensionlcss, time r = uj p t 1 co-moving position £ = u> p (t — z/v g ), longitudinal 
velocity = v z /c, Langmuir phase velocity p = v p /c = v g /c, vector potential a = a± = -^^A, 
electrostatic potential <f> = and number density perturbation p = (n e — tiq) /no, the continuity 

equation (JJJ can be written in scaled variables as 

J^ + ^[(i-/VA0(i+p)] = o, (ii) 

and Poisson's equation JHJ becomes 



^ = PpP- (12) 

Assuming an initially quiescent plasma in the initial absence of laser fields (i.e., rj, <f>, /3, a — * as 
r — ► — oo for fixed £ or as £ — > — oo for fixed r), the conservation of transverse canonical momentum 
implies that 

p_L = mca± — mca. (13) 
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The relativistic factor 7 then factors into transverse and longitudinal contributions as 

1 



7 = 7J-7|| = V 1 + \ a \ 2 



(14) 



Using these relations, we eliminate any explicit appearance of v± or p± in the longitudinal component 
of the momentum equation |JSj, thereby obtaining 



1 d 



1 1 d 



±m + (i - m) = - Jp oi^ - 2 ^\-\ 



(15) 



Now, we further make the Quasi-Static Approximation (QSA) (see, e.g., |Tl| and wherein we 

assume p = p(£) and (3 = /?(£), i.e., that the plasma response is independent of time r in the co-moving 
frame, implying that the plasma wave itself moves without dispersion at the fixed average group velocity 
P p of the driving lasers. The QSA also specifically requires that any distortion or change in the shape of 
the laser envelope remains negligible during the typical interaction time with any one transverse slice of 
plasma, which will be 0(T). This constraint is satisfied automatically in our model, since we have already 
assumed that laser envelope evolution is negligible throughout the entire interaction, which is of duration 
0(T + L/c), where L is the total plasma length. In the QSA, all r derivatives in equations describing the 
plasma response can then be neglected, and the continuity equation (|llf) can be immediately integrated, 
yielding, for our quiescent initial conditions, 



(l-/3/&)(l + p) = l. 



(16) 



Expressing \a\ 2 in terms of 7 and /3, the longitudinal momentum equation (|15f) can also be integrated 
to obtain: 



7 (i-/y?) -0 = 1. 

After some algebra, the Poisson equation l|12f> then becomes 



d 2 

^ = Phi 



Pp 1- 



7 P 2 (1 + ' 



(17) 



(18) 



This single, second-order equation describing the nonlinear plasma response is valid, within the QSA, 
for arbitrary relativistic electron velocities provided the Langmuir wave amplitude remains below the 
nonlinear wave-breaking limit determined by: 



> -1 



1p 



12 - 



-i + 2i. 

7 P 



(19) 



In the high phase-velocity limit appropriate to highly underdense plasmas, where < 1 — (3 P <C 1 and 
7 P 3> 1, the dynamical equation l|18|) may be Taylor expanded and further simplified to 



(20) 



Because we have effectively taken j p —> 00, this equation remains mathematically well-defined for any 
(p > — 1 and any \a\ 2 > 0, but the actual bound (|19|l remains a more physically trustworthy limit for 
wave-breaking. After solving for </>(£), we can determine the electric field using 



d_ 



(j) = j3 p £\\ 



(21) 



where = £\\{C) = E z /Eq is the scaled longitudinal electric field within the QSA. The density pertur- 
bation p(£) may then be determined using the scaled Poisson's equation i|12[l . the velocity /3(£) using 
the first integral (|16fl . and finally 7^) by using l|14f> . 



Under our assumptions, the normalized vector potential 



a(?,r) 



mc- 



■A (z=(t- £)v g /w p , t = r/ujp) 



(22) 
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is taken to be a prescribed function of time and position, describing the unperturbed laser fields of 
duration T throughout the plasma of length L, with a modulated envelope which is slowly chirped but 
travels along at a fixed group velocity v g . For otherwise arbitrary laser fields, \a\ 2 will still have both 
fast (carrier frequency) and slow (beat frequency) dependence on both r and £, although considerable 
simplification is possible for the assumed PBWA form, namely: two weakly detuned, slowly-chirped, co- 
propagating, nearly-plane- wave, flat-topped lasers with positive helicity. The normalized vector potential 
a, representing the beating drive lasers may then be explicitly written as 

a = ai + a 2 = \ [e+a^ 1 + c.c] + \ [e+aae 4 ^ 2 + c.c.] , (23) 

where e + = -j= (x + iy) . 

At the leading edge (z = 0) of the plasma, the laser phases tpj (j = 1,2) are given, for t > 0, by 



ipj(z 



Z 

= 0,t) = ^ jo - J dt'uijit'), 



(24) 



where the ipj are real constants depending on initial laser conditions, and the instantaneous carrier 
frequencies u>j(t) = Wj{z = 0,t) = — -^ipj(z = 0,t) allow for slow chirping of one or both of the lasers, 
such that 



i d . 



(25) 



and 



l w l — I < \u>2 — 0J2 I < \0Jl 



- UJ 2 \ 



\U)l - U>2\ 



10 



Pi 



(26) 



where we define the average carrier frequencies 

T 

1 

to 



I J dt' ujj (f) = (z = 0, t = 0) - ^ (z = 0,t = T)] /T, 





and the overall average carrier frequency is then defined as 

u> = \ (a>i + w 2 ) , 



(27) 



(28) 



Although taken as piecewise constant here, the amplitudes aj, j = 1, 2 could more generally also include 
suitably slow (i.e., slower than the plasma frequency, but possibly comparable to the chirp rate) time 
and position dependence as well, to model more realistic ramping of the drive fields and some 2D effects. 

Within the plasma (0 < z < L), we assume that the local laser frequencies uij(z,t) = —■g^'ipj and local 
wavenumbers kj(x,t) = -§^^pj > each satisfy the ID electromagnetic dispersion relation 



U3„ 



2 P i 2 7 2 

u = — + c k , 

7o 



(29) 



where the constant 70 > 1 parameterizes an average global nonlinear shift in the effective plasma 
frequency due to transverse electron quiver. For electromagnetic waves satisfying this dispersion relation, 
the group velocity v g is given by 



1/2 



(30) 



which in our undcrdense case may be approximated as 



V B (W) 



1 ^ 



2 7o^ 2 



(31) 
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We take as our reference group velocity v g , the expression (|30() evaluated at the average carrier frequency 



Vg = C(3 p = Vg(0j) W C 



w 2 - 



k(u> 2 ) — k(u>i) ' 



(32) 



which represents the characteristic velocity at which both laser envelope modulations and Langmuir 
phase-fronts travel. 

Because the relevant time-scales satisfy ujJ 1 -C w^ 1 <C T, the fast carrier oscillations in \a\ 2 at the 
harmonics 2 u>\ and 2 cj2 and the sum frequency u)% + lo 2 will average away, leaving only the slowly- 
varying contribution to \a\ 2 : 



a 2 = (\a\ 2 ) = \ \\a L \ 2 + \a 2 \ 2 + afae^-^ + ai a* 2 e~ 1 ^-^ 



(33) 



Since the group- velocity dispersion effects {-^i^> and higher-order terms) remain small in the underdense 
regime, linear propagation into the plasma results in a beat phase given by 



t-Z/Vg 

Tp 2 (z,t)-i> 1 {z,t) = Atlj + I dt' M*')-wi(t')] +0 



1 W p LOpL 

7o c 



(34) 



where Aipo is just a constant, equal to the difference of the initial laser phases. The neglected terms 
limit the validity of the constant group velocity approximation to interactions lengths L less than the 
so-called dispersion length L diap : 



uj 3 c 

Li S: J-i disp ~ — 5 . 

p V 



(35) 



For accelerator applications, the useful interaction length is already limited by the dephasing length L d 
beyond which accelerated electrons cannot gain energy from the electrostatic field: 



Up Up w 



Laisp *C L diap 



(36) 



so that our constant group velocity approximation imposes no further restriction on the interaction 
length. 

By appropriate choice of the initial phases, we may take both a\ and a 2 to be real and nonnegative. 
Defining the co-moving normalized beat frequency 



Aw(£) = [uj 2 {t = f/wp) -wx(t = €/wp)] /uj p , 



(37) 



the beat phase 



m = m + / df mo, 



(38) 



for ip(0) = AipQ, the normalized beat amplitude 



e = 0,10,2, 



the average normalized intensity per laser 



-2 1 r 2 1 21 
a = ± [af + a 2 \ , 



(39) 
(40) 



and the electron quiver factor 



7o 



(41) 
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the reference group velocity v g is fully determined, and the slow part of the normalized ponderomotive 
drive may be written as a function of £ only: 



a 2 (£) « (0 6 (u p T - [a 2 + ecos </>(£)] , 



(42) 



where 0(£) is the usual Heaviside step function. Using this form for the ponderomotive drive, the 
equation of motion (|2U|) may be written as an ordinary differential equation in the co-moving coordinate 

= 2 



1 + a + ecosV'(^) 

Jx + W 



- 1 



(43) 



with the initial conditions 0(£ = 0) = </>'(£ = 0) = 0, valid for r > and r — ui p L/v g < £ < r, while 
otherwise </>(£, t) = 0. Equation (|43|) is used for all numerical simulations discussed subsequently, and is 
the starting point for our analysis of autoresonance. 



III. HAMILTONIAN FORMALISM 



To study autoresonance, we now develop the Hamiltonian formulation of (|43(l . Our goal is an expression 
in terms of canonical action- angle variables, for which the phase-locking phenomenon is most readily 
analyzed. First, note that the dynamical equation for the electrostatic potential l)43|) can be derived 
from the Hamiltonian 

i a 2 + ecos^(Q . a 2 + e cos ?/■(£) , . 

+ 2(l + fl = U ^ P) + 2(1 + 0) ' (44) 



1 2 1 



n(<P,p;Z) = 1 2 



with the scalar potential regarded as the generalized coordinate, p = <p' = -^4> — j3 p E\\ regarded as 
the canonical momentum conjugate to 0, and £ taken as the time-like evolution variable. In this way, 
the plasma-wave dynamics are seen to be analogous to those of a one-dimensional forced nonlinear os- 
cillator. The component Hq of (|44|) represents the Hamiltonian of the free oscillator, involving one term 
T(p) = \p 2 analogous to the kinetic energy of the oscillator (which is in fact proportional to the electro- 
static potential energy density of the plasma wave), and another term corresponding to an anharmonic 
effective potential V{<j)) = \ [1/(1 + 0) + <f> - 1] . The remaining driving term AH. = H.(<f),p; £) — 7io(</>,p) 
corresponds to the time-dependent forcing of the oscillator. This forcing is external in that it depends 
on a prescribed function of £, but because of the nature of the relativistic nonlinearity, the effective 
strength of this forcing depends on both the intensity of the driving lasers and the instantaneous value 
of the dynamical variable <p associated with the excitation. Because of the ponderomotive nature of the 
coupling, the effective forcing is always positive, i.e., includes a constant as well as a purely oscillatory 
part. Both of these features differ somewhat from previously-studied autoresonant systems, and will 
have important implications for the dynamics. 

In the absence of forcing (i.e., e = a 2 = 0), the dynamics governed by 7io(<p,p) = T{p) + V(</>) are 
conservative, so the oscillator energy, i.e., the value H of the Hamiltonian TLq along any particular 
unperturbed orbit, remains constant. In the physically allowed region <j) > — 1, the effective potential 
V(4>) is nonnegative and possesses a single minimum V = at <f> = 0, while V(4>) — > +oo as cf> — > — 1 + 
or <j) — * +oo. So for any value of energy H satisfying < H < oo, there exists a phase space trajectory 
(0(£; H), </>'(£; H)j which is a closed periodic orbit, mirror-symmetric about its turning points, for which 
— 1 < 4>(£) < oo throughout, and which is unique up to its overall phase, or initial position at £ — 0. 

Making a canonical transformation to the action-angle variables of the free oscillator, <fi = <j>(T,8), 
p = p(T, 8), we can express (|44|l as 

where the action X is defined in terms the area in phase space contained within the unperturbed orbit 
H)^'(£-H)) of energy H: 

<!>+ 

1 =— I pd(j)= - /V #. (46) 



2ir T ir 
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Here, </> + and <f>— are, respectively, the upper and lower turning points of the orbit, at which V(</>±) = H 



b±=H± ^H 2 + 2H, 



(47) 



and in ((46(1 we have used symmetry to reduce the integration path to the segment where p > 0. By 
making the change of variables 4> — 4>+ ~ (4>+ ~ 4>-) sin 2 (u), the action 146(1 can be calculated with the 
help of a standard integral table, e.g., equation 2.584-13 on p. 163 of |38j : 



I = 



M 2 



tt/2 



V 1 + fa 



du 



sin 2 (it) cos 2 (it) 



3tt 



l+H- \JH 2 + IE 



1+0+ 
1/2 



■ sin 2 (w) 

(1 + H)E(k) 



l + H- VIP 



2H 



-K(k) 



(48) 



Here, K(k), E(k) are complete elliptic integrals of the first and second kind, respectively, whose modulus 
satisfies 



2(1 + H)\jH 2 + 2H - 2H(2 + H) 



1/2 



(49) 



At this point, one could in principle use 1(48(1 to find Hq(T), but fortunately such a cumbersome inversion 
will not be necessary. The normalized (i.e., dimensionless) nonlinear frequency Q,(H) of the unforced 
oscillator is given by 



dH _ (diy 1 _ tt [1 + H - VH 2 + 2H 
dl ~ IdHJ ~ 2 



1/2 



E{k) 



(50) 



To put 1(44(1 in the desired form 1(45(1 . the remaining ingredient needed is the canonical transformation 
<j> = 4>{I, 0) . An explicit formulation of this will also not be needed, and we may proceed by formally 
assuming that we have made this substitution. Then, 4> = <j>(I, 9) is a periodic function of 6, and we can 
expand the driving term appearing in 1 (45(1 in a Fourier series as 



:cosi/'(£) 



[a 2 +e cos ?/;(£)] ]T b n (X)t 



(51) 



2 [1 + ^(1,6)] 

Because ((51(1 is a real- valued function, the Fourier coefficients must satisfy b n — b*_ n , and are defined by 



b n {T) 



d9- 



2irJ ' 2[1 + </>(!, 6>)]' 
o 



(52) 



We put ((52(1 into a form more amenable to calculation by changing variables using 9 = il(H)£, and 
integrating over one period < £ < H of the co-moving coordinate, where S = S(iJ) = 2n/D,(H). 
Furthermore, if we choose the origin of the orbit associated with the energy H, or equivalently with the 
free action 1(H), such that (/>(£ = 0;H) — </> + , the potential <fi(£;H), and hence [1 + </>(£; H)] -1 , are 
symmetric about the point £ = 3/2. Since the imaginary parts of ((52(1 are obtained by integrating over 
the antisymmetric functions sin(nf2£), the 6„'s are purely real and given by: 



E(H) 



b n {H) = b n (H)* = b- n (H) = b- n (H)* 



E(H) 



r/£ 



cos [nn(H)£\ 
2[l + 0(e;fl)]" 



(53) 



Up to this point, no approximations have been made in the canonical transformations to the action-angle 
variables of the unforced oscillator. Now, we make the Single Resonance Approximation (SRA) [33 to 
1451) , by assuming that the rapidly oscillating terms of the Hamiltonian average to zero and contribute 
negligibly to the dynamics. Anticipating the developments of Section llVI we realize that under certain 
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frequency-locking conditions derived there, autoresonant excitation occurs, wherein the plasma wave 
amplitude, or equivalently free-oscillator energy H(£) = X (</>'(£)) + V (</>(£)) , grows secularly such that 
the dynamical frequency f2(£) °f the forced oscillator approximately matches that of the unforced oscil- 
lator at that same amplitude, f2(£) ~ f2(ff(£)). Simultaneously, this instantaneous oscillator frequency 
follows that of the driving beat frequency of the wave, fi(£) ~ Au/(£). Under these assumptions, we can 
consistently neglect all terms in the sum Ij51|) except the constant term or those terms with frequency 
dependence ~ ± — Aw(^)] . The averaged Hamiltonian then becomes 

H(X, 6: = H (X) + eh(X) cos[9 - + a 2 b (X). (54) 



We now make an explicitly ^-dependent canonical transformation to the rotating action-angle variables 
(X, W). Using the mixed- variable generating function 

F 2 (±,0;O = [0-^(0]^ (55) 
our old and new coordinates are related by 

*=-# = »-^(0; (56a) 
oX 

X=^=X. (56b) 
Dropping carets from the new action, the Hamiltonian in the rotating frame is 

K(I,*;t)=K(I,9,i) + -Q£ = Ho(X) - A^(£)Z + eb x {X) cos * + a 2 b (l). 
The resulting SRA canonical equations of motion are 

9T = W)-^(0 + e-co S * + -a^ (57a) 
±X=-^ = eb 1 (X)sin% (57b) 
for which we next determine the required conditions for phase-locking. 



IV. AUTORESONANT RESPONSE 



The essential ingredients for autoresonance, as exploited in the present scheme, are: a nonlinear, oscilla- 
tory degree of freedom (in our case, the plasma wave) evolving, in the absence of any forcing, within an 
integrable region of phase space; a continuous functional relationship between the nonlinear frequency 
and energy of the oscillation, possessing a well-defined linear limit; an applied oscillatory driving force 
(in our case, the modulated ponderomotive envelope of the lasers) which is sufficiently small as to be 
considered perturbative, so that the notion of the nonlinear frequency of the unforced oscillator re- 
mains meaningful; and initial conditions and forcing profile consistent with Adiabatic Passage Through 
Resonance (APTR) - namely, an initially unexcited system (quiescent plasma), and an initial drive 
frequency sufficiently far from the linear resonance, with subsequent time-dependence that is sufficiently 
slowly- varying but otherwise arbitrary. 

The key consequence of the autoresonant beat-wave generation is the robust entrainment between the 
three relevant (normalized) frequencies: the beat frequency Aw(£) of the driving lasers, the instantaneous 
frequency f2(£) of the driven plasma wave, and the nonlinear frequency £l(H) of the unforced plasma 
wave. Assuming this phase-locking is achieved, amplitude control of the plasma wave can be simply 
understood via l|50() : because the frequency of the freely-evolving nonlinear plasma wave is a function 
of the energy, phase-locking of the driven wave to the envelope such that Aw(£) w ~ Sl(H) implies 
that changing the drive frequency will correspondingly change the oscillator energy. In our case, (I5UI) 
indicates that fi(if) is a decreasing function of the energy, so that in order to increase the plasma 
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wave amplitude one must decrease the beat frequency as a function of the co-moving position £ (or, 
cquivalently, in time t at the source of the driving lasers.) 

While we have indicated how autoresonance can lead to large plasma waves, we have not yet shown 
under what conditions such phase-locking occurs. To answer this question, we first consider the linear 
and weakly nonlinear response, valid up to the point where the normalized sweeping drive frequency is of 
the order of the normalized linear resonance, fl(H) « Aw(^) > 1. Next, we consider the fully nonlinear 
case, for which we derive adiabaticity requirements for autoresonance. 



A. Small Amplitude Response and Phase-Locking 



When the drive is first applied with its frequency above the linear resonance, the plasma wave amplitude 
(and, hence, H and 1{H)) are small and we can linearize equations l|57al - l57b|) . In this limit, the 
oscillator is harmonic, with Q(H) = 1, <fi(£) = \2H cos(£), and H = I, so that bo = —\H, b\ = \2H /4, 
and equations j|57J) become 

^ = i_Ia 2 -A^) + ^-^cos*, (58a) 

Note that the (1 — \a 2 ) contribution to corresponds, in normalized units, to the leading order 
expansion of the effective plasma frequency in the EM dispersion relation: 



u> Poit = ujp/jo = Up/y/l + a 2 s» u p (l - \a 2 + . . .) , (59) 

which is shifted from the bare value oj p due to the transverse quiver motion of the electrons in the 
applied laser fields. While we never explicitly invoked any small a 2 approximation, only the leading-order 
correction appears because, in making the SRA, we ignored terms of the form a 2 e m0 for \n\ > 1. But for 
sufficiently large intensity a 2 such terms can appreciably effect the motion despite being off-resonance. 

1/2 

One could, a posteriori, partially correct for this defect by replacing (1 — \a 2 ) with (l + a 2 ) in 



the equation of motion (|58a|l . but this will be unnecessary for the small values of a considered here. 
Physically, we should simply ensure that the drive frequency begins sufficiently far above, and then is 
slowly swept past, the effective (or renormalized) frequency w Poff , rather than the bare frequency lu p , 
if the difference is not negligible. Precise knowledge of the exact value of the effective linear plasma 
resonance including the effects of transverse quiver is not needed. 

For the driven plasma wave, we have Acj(£) ~ ~ 1, while we seek solutions for which |Au;(£) — 
< 1 as a result of special initial and forcing conditions: an initially unperturbed plasma, (/>(£ = 
0) = </>'(£ = 0) = 0; an initial tuning of the beat frequency above resonance, i.e., Aw(^ = 0) > 
tt>p off , and subsequently a slow downward frequency chirp through resonance, where the chirp rate is be 
characterized by a parameter a = a(£) = — ^Aw(£), with < a -C 1. 

For a linear frequency chirp around the effective plasma frequency, 

Aco(£) = 1 - \a 2 - at£, (60) 

the simple harmonic oscillator equations (|58l) have analytic solutions in terms of the Fresnel Sine and 
Cosine integrals |40j . We briefly summarize the extensive characterization of these solutions found in 
|3lj . When the drive is first applied far from resonance, the oscillator response can be divided into two 
components: one ringing component precisely at w Poff , and the other at the driving frequency Aw(0), 
both of small amplitude. The singular term ~ I -1 / 2 in (|58a|l allows for a large change in phase at small 
amplitude without violating the requirement that </>(£) remain smooth, so that the response at the driven 
frequency can adjust itself to the drive, and phase-locking can occur. As the frequency is swept toward 
the resonance, these driven, phase-locked oscillations grow. Meanwhile, the response at the resonant 
frequency has no such phase relation with the drive, and remains small. In this way, we essentially have 
one growing, phase-locked plasma wave when the drive frequency reaches the resonant frequency. 
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FIG. 1: Demonstration of the critical autoresonant behavior for e = 0.005. In (a), chirp rates 1% below critical 
(solid line) yield final longitudinal fields above linear wavebreaking, whereas chirp rates 1% above saturate at 
0.6-Eq- (b) compares the phase locking in autoresonance (solid line) to that of non-autoresonant behavior (dashed 
line). The critical drive (mixed line) is included for reference. 



As the resonance is approached, the amplitude of the plasma wave begins to become large and one 
must account for the growing nonlinearities. We therefore expand the expression for the free nonlinear 
frequency (| 5 f )|> to first order: 

n(H) = 1 - ftf = 1 - fx, (61) 

and continue to use the linearized frequency chirp i|60|) . Making the change of variable A = 4-V/2X, we 
have 



d 



— \$> = a £ - 2§6^4 2 + — cos*, (62a) 



-^:„4 = esin*. (62b) 

This set of equations can be reduced to a single first-order ordinary differential equation by defining the 
complex dynamical variable Z = —^256/ 3a.Ae l * , a re-scaled independent variable £ = ^/a^, and the 
dimensionless parameter /i = e-y/3/(256o! 3 / 2 ), to obtain 

i-^Z + (C - \Z\ 2 ) Z = [i. (63) 

Thus, the weakly nonlinear problem is now described by a dynamical equation with a single parameter 
fx that combines the drive strength e and the chirp rate a. It has been found numerically |41| that the 
solution to (|63|l has a bifurcation at the critical value fj, = fi c « 0.411. For /i < /i c , the plasma wave 
response quickly dephases from the drive, resulting in only small excitations. In contrast, for /i > /i c , 
phase-locking occurs and the plasma wave can grow to large amplitude. This critical behavior with 
respect to fi translates into a critical drive strength e and chirp rate a for the nonlinear oscillator to be 
autoresonantly excited. The condition is 

Thus, for a given laser intensity one can readily find the maximum chirp rate that can be tolerated and 
still obtain high amplitude plasma waves. We demonstrate the sensitivity of this critical behavior for 
e = 0.005 in Fig.^ which plots numerical solutions to the full quasi-static equation of motion l)43|l . Plot 
(a) shows the excited longitudinal electric field for chirp rates just above and just below the critical rate 
a c « 1.28 x 10 -4 ; Plot (b) demonstrates the dynamic frequency-locking that occurs in autoresonance, 
and compares this to the case where autoresonance fails to occur. 



B. Fully Nonlinear Autoresonant Response 



If phase-locking is maintained through the weakly nonlinear regime, the amplitude continues to grow 
and one must consider further nonlinearities beyond those included in l|63|) . In this case, there arises 



14 



more stringent restrictions on the chirp rate for adiabatic phase-locking to persist. We calculate this 
condition by first finding a second-order equation for the phase ty: 



d 2 



d_d® _ d 2 JC dJC 

d£ d£ ~ dld^ dl dl 2 dl ' dld£ 



d 2 JC dK, d 2 JC 
+ 



d_M 

d£ dC 



(65) 



where we have made use of the usual canonical Poisson bracket with respect to the rotating-framc 
action-angle variables: 



{F,G} 



d£dG 
9* dl 



d£dG 
dl <9*' 



(66) 



Using this expression ffity as well as l|57a|) . the phase \& is seen to obey the following second-order 
equation: 



d 2 db d 

_v£r + e -± sin W— W + Au; - eh (1) sin * 



an 

dl 



_ 2 d 2 b , 9 2 6i 



dl 2 



+ e 



dl 2 



cos 



= 0. 



Now, we assume (see, e.g., (43) that the free action can be written as 

1 = In + AX, 



(67) 



(68) 



where 1q = 2o(£) is the slowly- varying, secularly-evolving action about which there are small oscillations 
given by Al = AZ(£). These oscillations correspond to fluctuations in \& about its (slowly- varying) 
phase-locked value ^> = *(^), an example of which can be seen in Fig. db). In the autoresonant case 
(solid line), we see that as the plasma wave is excited, its frequency does indeed make small oscillations 
about the drive frequency. We further note that as one decreases the chirp rate a from its critical value 
a Cl we obtain more total oscillations in frequency over the longer excitation time, but with a slightly 
smaller magnitude of excursions from the drive frequency. 

Using the form (|B^|l for the action, the lowest-order equation for the phase is identical to (|BTjl . with 
1 being replaced everywhere by In. In this way, the phase itself is seen to obey a nonlinear oscillator 
equation, with an effective nonlinear damping (or anti-damping) term, and a conservative "forcing," 
described by the terms in the braces, which is derivable from an effective "potential" (not to be confused 
with any previously-mentioned potential) whose shape is dictated by the slowly-evolving action 1 and 
the drive parameters e and a. Phase-locking then corresponds to trapping of $ in a basin of this 
effective potential. In order for the phase \E' to be trapped about some value the effective potential 
must possess a local minimum there, and the non-conservative term must either provide damping or 
else remain sufficiently small if excitatory. In fact, the non-conservative term is expected to be small 

compared to the second term in the conservative force. Their ratio is given by (j^f^/J^) (^f) ' 



l g6i 

hi dSl 



1_ 

<> 



O(l), while ^ <§; O(l) because the nonlinear phase 



where typically (^/^) 
oscillations are slow compared to u p . Actually, over most of the typical range of parameter values, the 
last two terms in the conservative forcing are small compared to the first two terms, since they are 
higher order in the drive strength, and as a first approximation the phase oscillations are governed by 
the "biased" pendulum equation 



, . dn 

e6i(Xo)^sin*, 

Oln 



(69) 



although we will continue to work with the full equation (|67|l . Clearly, for a given 1 , if the normalized 
chirp rate a remains sufficiently small compared to the normalized drive strength e, then the effective 
potential will be of the tilted-washboard variety, with a series of periodically-spaced local minima in ^ 
at intervals of 2tt. As a increases, the depth of these wells decreases, until they finally disappear, as does 
any opportunity for phase-locking. 

Thus, a necessary condition for trapping is that the effective force can actually vanish at some fixed 
point 'J: 



<*(£) +e6i(To)sin* 



dn 

92o 



>d 2 b 



d 2 h - 

+ £ COS V 



dl 2 dl 2 



= 0. 



(70) 
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FIG. 2: Shows the maximum plasma wave amplitude obtainable before the "slowness" condition 1711 is violated. 
In (a), we plot the maximum longitudinal field as a function of the chirp rate for four different driving strengths. 
Dotted lines indicate hysteresis at a c given by 16411 . (b) compares the theory (line) with numerically determined 
saturation for e — 0.005 



For given To, e, and a, this determines the average, slowly-varying phase W about which trapping occurs. 
In the linear (In — > 0), weakly forced (0 < e <^ 1) or weakly chirped (a — * + ) regimes, this phase 
value is known [31| to be ^ ~ ir, but for nonlinear plasma waves, stronger forcing, or faster chirping, 
this phase can shift appreciably. Since |sin^|, |cos^| < 1, equation 1)70(1 also imposes an upper bound 
on the frequency chirp a(£) for which the phase can remain trapped in autoresonance, regardless of the 
actual value of the phase at which locking occurs. Setting |sin^| = |cos'F| = 1 above, we obtain an 
upper bound on a beyond which any phase-locking is impossible: 

< a(0 < e|&i(X )| 

Again, for realistic parameters, the force balance typically resides predominately between the first two 
terms in (70}, and hence this upper bound, although approximate, is expected to provide a reasonably 
tight cutoff for autoresonant phase-locking, which has been confirmed by numerical simulation. This 
inequality can also be thought of as giving the maximum achievable plasma wave amplitude (implicitly 
as a function of In) for a given chirp rate and laser power. We show the dependence of the saturated 
longitudinal field on the chirp rate as a solid line in Fig. ffia) for a number of different drive strengths. 
For a fixed drive strength e, the maximum attainable electric field jumps discontinuously at the critical 
chirp rate a c given by l|64|). The dotted lines correspond to solutions of (|71ll that cannot be accessed 
when starting from vanishing initial longitudinal field, and as such constitutes a form of hysteresis in 
the excitation. The solid lines show the stable branches for the case of interest: excitation from a 
quiescent plasma via slow downward chirping from a drive initially above resonance. Fig. Gib) shows 
a comparison of the theoretical maximum amplitude and that found by numerically integrating the 
quasi-static equation of motion l|43|l . 

Asymptotic expansions and numerical plots for 0,(1) and b n (T) reveal that the right-hand-side of H71|) 
actually decreases with increasing In (at least beyond moderate values of In), making the adiabatic 
nonlinear phase-locking condition l|71l) become increasingly stringent as the plasma wave grows. For 
fixed chirp rate and drive amplitude, the growing plasma wave will eventually fall out of frequency- 
locking. In the physical pendulum limit discussed above, note that when this cutoff is reached and 
autoresonance is lost, the average phase should be near ±7r/2. Mathematically, this saturation can be 
traced back to the nonlinear nature of the forcing in Q43[l. where the effective strength of the forcing 
depends both on the external drive intensity and the plasma wave amplitude; physically, to the fact that 
as the plasma wave grows, the electrons become increasingly relativistic and therefore less susceptible to 
displacement via the ponderomotive forcing at a fixed intensity, so the effective strength of a given drive 
actually decreases. In principle, this may be counteracted by either slowly increasing the laser amplitude 
or, more practically, by slowly decreasing the chirp rate in time. In the latter case, the rate of growth of 
the plasma wave will fall off because it is directly tied to the decreasing drive frequency, but the wave 
can remain frequency- locked with the drive indefinitely (with some oscillating excursions.) 

As a result of relativistic detuning, recall that in the original RL/TD scheme, the plasma wave amplitude 
exhibits slow (compared to Wp 1 ) nonlinear modulations which appear as beating, i.e., periodic amplitude 
oscillations up to the RL limit and back to a nearly unexcited state. As the wave is driven to a high 



on 


+ a 2 










dl 




dll 


+ e 


dT 2 





(71) 



1G 



amplitude its phase slowly slips until it is more than ~ ir/2 out of phase with respect to the laser beat and 
then gives its energy back to the lasers, then continues to shift further out of phase, only to be re-excited 
when its amplitude approaches zero and it can re-establish phase matching with the drive. Both the 
frequency difference between plasma wave and drive and the phase-lag exhibit continuous oscillations in 
time. 

In the DMG scheme, which relies on autoresonance but starts at linear resonance, the plasma wave 
amplitude not only can peak at a higher maximum than in the original scheme, but typically will sustain 
a higher average value at long times, exhibiting a nonlinear ringing about some non-zero saturated value 
rather than a full beating. We suspect that this behavior results from two features of the chirping: when 
the autoresonant adiabatic condition fails, the wave phase is closer to its neutral value, with respect to 
energy exchange with the drive, and after autoresonance is lost, the frequency difference and phase lag 
each grow secularly in time rather than exhibiting oscillations. Depending on initial parameters, the 
initial growth up to the absolute maximum can either be essentially monotonic, or exhibit a "staircase" 
behavior, with intermittent plateaus or dips at local maxima between periods of resumed growth, before 
finally leveling off with some ringing. When starting from resonance, the originally published numerical 
examples 27j and our own simulations suggest that the observed ripples in the amplitude excitation are 
actually minimized at some intermediate value of the chirp rate, implying that the excitation may not 
be fully autoresonant in some sense. 

In our autoresonant scheme based on APTR, the behavior more closely resembles that in DMG scheme, 
but exhibits more nearly monotonic growth, higher peak fields, and less ringing after saturation. Specif- 
ically, some slow oscillations in the amplitude may be present early on, but once nonlinear phase-locking 
has been achieved, the growth remains monotonic or virtually so until finally the amplitude levels off 
and appears to almost saturate, with only a very small amount of subsequent ringing. 

The extent of the amplitude and phase excursions, which affect how regularly the excitation grows and 
how closely the phase is locked during autoresonant excitation, and also how much ringing persists after 
saturation, are determined by how deeply the phase is trapped in its effective potential well. This in 
turn depends both on how steep and how deep is the available well (determined by the plasma wave 
amplitude and drive lasers parameters), and to what extent the phase can be nudged into position near 
the bottom of the well and kept there (determined by the initial conditions and the adiabaticity of the 
chirped forcing). 

Numerical simulations suggest that, as one would expect, the extent or depth of this trapping is improved 
by using a stronger drive (at least up to some moderately strong value), starting the drive frequency 
further above resonance, and chirping more slowly. In practice, of course, each of these strategies 
involves trade-offs. Increasing the drive strength increases the growth rates for laser-plasma instabilities 
that might disrupt the forcing. Either increasing the initial frequency up-shift or decreasing the chirp 
rate decreases the final amplitude that can be reached during a fixed interaction time. 



V. EXPERIMENTAL CONSIDERATIONS 



Unfortunately, as has been alluded to previously, the PBWA does not have unlimited time to be excited, 
as dclctcrous instabilities will eventually destroy wave coherence. For the parameters of interest, the 
oscillating two-stream (also referred to as modulational) instability limits the lifetime of coherent Lang- 
muir waves to the ion time-scale, i.e., for times of the order of a few l/wf. Although it is possible that 
the growth of this instability may be mitigated somewhat by the use of a chirped laser, in this paper we 
will use as a conservative figure the results of Mora et. al. [T^ | to set the time limit during which we 
can excite a coherent plasma wave suitable for accelerator applications. 

For the relatively cold plasmas and moderately intense lasers we consider, it is shown in |l3j that the 
growth rate of the oscillating two-stream instability is approximately equal to Wj, and that this instability 
impedes plasma wave excitation and destroys coherence after about 5 e-foldings. Thus, we see that the 
drive lasers should have time duration T < 5/a>j. If one chirps the drive frequency leading to a total shift 
6u> during the autoresonant excitation, then the normalized chirp rate is limited to 

, . . ui„ Suj . . . 
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FIG. 3: Plasma wave excitation for a 10/im CO2 laser with intensity 2.7 x 10 14 W/cm 2 (e = 0.02). In (a), the 
APTR case has to(0) = 1.15oj p , and = uj p at i ~ 40ps (a — 0.00065). Total chirp is 1.5% of laser frequency. 
For comparison, we include the envelopes with no chirp, and DMG chirp starting on resonance, (b) is the same 
as (a), with u> p changed by 10%. In (c), we show uniform accelerating fields via APTR for variations in uj p of 
order 10%. (d) demonstrate robust field excitations of 5-10GeV/m for density errors of order ±20%. 



or approximately a > 2.3 x 10~ 3 5oj/u) p for singly-ionized Helium. Below, we choose two experimentally 
relevant parameter sets, one corresponding to a 10 /jm CO2 laser; the other, to a 800 nm Chirped Pulse 
Amplification (CPA) [43l |44| TkSapphire laser system. We demonstrate how, beginning with the laser 
frequency above the linear resonance and then slowly decreasing it, one can robustly excite plasma waves 
to amplitudes larger than the cold, linear wave-breaking limit in times commensurate with onset of the 
oscillating two-stream instability. 



A. CO2 Laser at 10 /im 



We consider parameters roughly corresponding to the most recently published UCLA upgrade |l8j . 
We assume two pulses of duration T = 100 ps which enter the plasma at t = 0, central wavelengths 
near A = 10 /im, and normalized intensities a\ — 0,2 — 0.14, corresponding to a normalized drive 
strength e = 0.02, so the threshold condition 164[) implies that the normalized chirp rate should satisfy 
a{t) = — -^Aui(t)/uj p < 0.0009. We choose a linear chirp, so that in physical units the beat frequency 
is given by Auj(t) = (fi a + auj p t), where fio — 1-15 and a = 0.00065, with a total frequency sweep from 
Aoj(t = 0) = 1.15 to Au>(t = T) = 0.74. For these parameters, we collect the results from simulations 
integrating the quasi-static equation of motion (|43|l in Fig. |3J Ela) demonstrates the excitation of 
a uniform accelerating field E z of 10 GV/m, which is above the linear cold wave-breaking limit of 
Eq ~ 8.8 GV/m, but below the cold relativistic limit £ WB ~ 61 GV/m. The total chirp is modest, only 
about 1.5% of the laser carrier frequency 2ttc/X. 

For comparison, we also plot the simulated envelopes of the longitudinal field for the resonant RL/TD 
case Aui(t) = Aw(0) = 1, and for the chirped DMG scheme starting on linear resonance, Auj(t) = 
(1 — aujpt) , but using the same chirp rate as above. The resonant case demonstrates the characteristic 
RL limit of E z < E RL = (16e/3) 1 / 3 £ ; o ~ 4.2 GV/m, whereas the DMG scheme fails to achieve appreciable 
dynamic phase-locking, and the final plasma wave amplitude is about the same as in the resonant 
(unchirped) case. Using approximately these parameters, UCLA experiments have inferred accelerating 
amplitudes up to 2.8 GV/m over short regions of plasma. More recently, plasma density variations 
corresponding to E z w 0.2 — 0.4 GV/m have been directly measured with Thomson scattering |45| . 
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FIG. 4: Plasma wave excitation for a 800 nm Ti:sapphire laser with u/lo v = 25, intensity 2 x 10 17 W/cm 2 
(e = 0.09). APTR parameters are Aoj(0) = 1.2, and Alu = 1 at t w 1.4 ps (a = 0.00225). Total chirp is 3% of 
the laser frequency. In (a), we see that the DMG and APTR scheme give approximately the same final fields, 
about three times that for on-resonance. In (b), we changed u v by 10%, and the excitation for the DMG and the 
on-resonant scheme drop considerably, while APTR maintains its large excitation, (c) demonstrates that APTR 
gives robust, uniform accelerating fields for wide variations in density, while (d) indicates the insensitivity of 
APTR to errors of ±35% in the plasma density, while still yielding large amplitude plasma waves. 



Perhaps more important than the higher amplitude field in the autoresonant APTR case, is the fact 
that excitation is very robust with respect to mismatches between the beat frequency and the plasma 
frequency. In practice, these mismatches inevitably result from limited diagnostic accuracy or shot-to- 
shot jitter in the plasma or laser parameters. Because one sweeps over a reasonably broad frequency 
range and one only needs to pass through the resonance at some indeterminate point during the chirp 
history, no precise matching is required, and the exact value of the plasma density need not be accurately 
known. This robustness is demonstrated in Fig. Etb)-(d). Plots (b) and (c) show the longitudinal field 
profile attained when there are variations in the density, and we see that APTR yields uniform, large 
amplitude fields over a wide range of densities. In Fig. E^d), we plot the final accelerating gradient 
achieved via APTR when we vary the value of ui p over a range of ±10%, from its "design" value, 
while keeping the laser parameters fixed. We see large levels of excitation for a wide range in plasma 
variation, corresponding roughly to density mismatch/errors up to 20%. Thus, not only is autoresonant 
plasma wave excitation effective in avoiding saturation from detuning, it also mitigates experimental 
uncertainties in or shot-to-shot variations of plasma density. 



B. Ti: Sapphire Laser at 800 nm 



Here, we analyze a representative case for a TkSapphire CPA laser in a singly-ionized He plasma, with 
no = 1-4 x 10 18 cm -3 , so that ui p /ui = 1/25, and a laser duration T = 3.2 ps, chosen to correspond 
to the modulational instability limit. If we consider two 1 J pulses compressed to this time duration 
and focused to a waist of wo = 6 /im, this implies intensities of Io = 2.0 x 10 17 W/cm 2 , so that, with 
A » 800 nm, we have a x = a 2 = 0.3, and e = 0.09. We choose Auj{t = 0) = 1.2, Aw(t = T) = 0.5, so 
that a — 0.0025. The resulting plasma wave excitation is shown in Fig. 0fa). Here, we see maximum 
longitudinal electric fields E z sa 260 GV/m, corresponding to ~ 1.6 Eo w 0.25i? WB - For comparison, we 
also plot the resonant case, for which detuning results in maximum fields corresponding to the familiar 
RL limit (lee/S) 1 / 3 ^ « 125 GV/m, and the DMG case (starting on resonance), which yields results 
similar to the APTR case. The distinction between passing though resonance and starting on resonance 
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can be seen, however, in Fig. H£b), which indicates that a small change in the value of the plasma 
frequency has a dramatic effect when starting on resonance, but little effect when passing through 
resonance. The uniform, robust acceleration fields obtained via APTR are shown in Fig.Efc) for density 
variations ±20%. Finally, Fig. 0^d) shows the robustness of autoresonant excitation, for which density 
imperfections of ±35% have little effect on the accelerating gradients achieved. 



VI. DISCUSSION: COMPARISONS, SCALINGS, AND EXTENSIONS 

In comparison to other PBWA schemes, including the fixed beat-frequency approach, either at (RL/TD) 
or below (TSS) linear resonance, the chirped (DMG) scheme, involving downward chirping from reso- 
nance, or the non-resonant PBWA, scheme, recently proposed by Filip et al. |45ll46|. involving strongly 
forced waves at frequency shifts well below resonance in a marginally underdense plasma, the autoreso- 
nant/ APTR PBWA enjoys a number of advantages, in terms of plasma wave amplitude, robustness, and 
quality. In previous sections we have seen how, for given drive laser intensity, autoresonant excitation 
yields longitudinal fields that can be considerably higher than the RL limit set by relativistic detuning of 
the plasma wave. We have also seen how APTR, i.e., slowly sweeping the frequency downward through 
resonance, provides a much greater degree of robustness to density mismatches, since neither the final 
amplitude nor frequency of the plasma wave is very sensitive to the precise location of the actual linear 
resonance, or to the precise chirp history. 

Direct comparison with the chirped DMG scheme is slightly more complicated, as revealed in Fig. 
because both schemes rely on autoresonance. In principle, given unlimited excitation time in the APTR 
case, and for any fixed drive laser intensity, the plasma wave can be autoresonantly excited to any 
amplitude at or below the nonlinear wave-breaking limit by choosing a sufficiently slow chirp rate. This 
is not always true for the DMG case starting at resonance, where performance appears to peak at some 
intermediate chirp rate, with rapid detuning for significantly faster chirp rates, and excessive ringing 
for significantly slower chirp. Without time constraints, the autoresonant APTR approach can always 
produce higher longitudinal fields for the same drive laser intensity or comparable longitudinal fields 
with smaller intensity. But in practice, the time allowed for excitation is inevitably limited, typically by 
ion instabilities as previously addressed, or even if these are somehow controlled, then by laser scattering 
or modulational instabilities, or ultimately by Landau or collisional damping of the Langmuir wave or 
hydrodynamic expansion of the plasma. 

If the DMG and the APTR schemes are compared for realistic parameters using the same drive laser 
intensity and chirp rate, then autoresonant APTR excitation consistently results in larger electric fields 
for excitation times on the order of few ui^ 1 . But if the DMG scheme is started precisely on resonance, 
but using a slower chirp rate so as to achieve roughly the same final drive frequency as the APTR 
approach (which started above resonance), then the DMG approach can experience phase- locking over 
a longer time and can achieve slightly higher fields in the finite excitation time. Essentially, by using 
a sufficiently strong drive, chirp rates that are just slow enough to be adiabatic, and an initial beat 
frequency precisely at resonance, the DMG trajectory can become autoresonantly phase-locked, but 
without having to waste time by chirping down from some point well above the resonance as is done in 
the APTR case. The catch is that such phase-locking behavior starting from resonance is quite sensitive 
to the initial conditions, and DMG will fail to achieve persistent phase-locking if the resonance is missed 
by just a few percent, so robust phase-locking behavior is unlikely to be experimentally reproducible, 
unless one begins with the beat frequency safely above resonance and adiabatically sweep through it, 
rather than trying to start precisely on resonance. 

In order to excite larger electric fields at a given plasma density, autoresonant/ APTR PBWA requires 
either more intense drive lasers or longer excitation times or both. Even if technologically feasible, moving 
to larger a 2 introduces its own problems (as discussed below), and in any event tends to undercut one 
of the primary comparative advantage of the PBWA over either the standard short-pulse or the self- 
modulated LWFA. Since excitation times in the PBWA tend to be limited by deleterious ion instabilities 
(oscillating two-stream instabilities and perhaps ion acoustic instabilities) growing on the w" 1 time-scale, 
while for fixed drive laser intensity both the maximum chirp rate and the wave-breaking field scale with 
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FIG. 5: Comparison of three plasma beat-wave schemes for the CO2 laser of Sec. IV Al Autoresonant/ APTR 
excitation is compared to the DMG scheme with identical chirp, and with a chirp chosen to maximize excita- 
tion. For identical chirp, DMG results in a field of approximately 50% less amplitude. The optimized DMG 
is comparable in amplitude to APTR, but a mismatch in plasma density of 2% has a sizable effect, unlike the 
autoresonant/ APTR robustness to density fluctuations of order 10% to 20%. 



the (electron-dominated) plasma frequency uj p , from the ratio 

we observe that, for given plasma density, the time available for plasma wave excitation may be increased 
simply by (singly) ionizing a gas comprised of a higher-mass atomic species. With heavier ions, one has 
more time to excite larger fields with the same laser power, or comparable fields with less laser power 
before ion motion disrupts coherence. In addition to ionic instabilities, laser instabilities are also of 
concern. Raman backscatter, with growth rate T RBS ~ \a^Junxi v and Raman forward scatter, with 
growth rate r RFS ~ a^f- , as well as related 2D self- focusing or modulational instabilities 0], are 
almost certainly not negligible for parameter regimes of interest, and their effects should be investigated 
in a more realistic simulations. There is some hope even that Raman scattering might be turned to 
our advantage: it may also be possible to exploit autoresonant excitation in counter-propagating or 
colliding-pulse schemes [47l l48l l49l | , where the enhanced growth of the plasma waves associated with 
Raman scatter may allow shorter excitation times for fields approaching the wave- breaking limit. 

Strictly speaking, our simplified dynamical model has demonstrated the robustness of autoresonant 
excitation only with respect to global density mismatches, and to small changes in the global chirp 
rate, drive intensity, or initial detuning. But we believe that some of this robustness should persist 
in the presence of moderate spatio-temporal variations and non-uniformity within the plasma and of 
complex spatio-temporal structure and dynamics in the lasers, features which are not captured by the 
idealized model used here. Effects such as ionization dynamics, thermal fluctuations, hydrodynamic 
expansion, electron and ion plasma instabilities, ponderomotive blow-out, and of course the nonlinear 
plasma excitation itself all lead to highly nontrivial plasma density profiles that vary in time and space 
and change the local conditions for resonance. Likewise, diffraction, tunneling ionization and the resulting 
ionization-induced depletion, refraction, and blue-shifting, as well as nonlinear focusing effects due to 
inhomogcncities and nonlinear back-action of the plasma wave on the EM fields, as well as Raman, 
self-modulation, and other instabilities, all can lead to appreciable distortion of the driving laser fields 
during the excitation process. The importance of these details for autoresonant wake excitation need to 
be further studied in the context of more detailed laser-plasma models. 

Any such variation, inhomogeneity, or complex dynamical structure in the evolving laser fields or plasma 
is likely to prove deleterious in the absence of autoresonance with adiabatic passage through resonance, 
drastically reducing wake excitation and uniformity. Longitudinal fields in some spatial regions may be 
resonantly excited to moderate or large amplitudes while other regions may be mismatched in density and 
experience very little excitation, or may be first excited and then subsequently de-excited as a result of 
phase slippage and energy exchange back to the lasers. As a result, we expect the final plasma wave to be 
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a highly irregular accelerating structure. But our autoresonant approach enjoys an intrinsic insensitivity 
and persistence, due to the local nature of the phase-locking. Provided only that the magnitude and 
scales for the non-uniformity are such as to allow an eikonal treatment of the waves, we expect [sij 
that at each position, the local plasma response will be autoresonantly excited by the drive, based on 
the local, slowly-varying values of the plasma frequency, drive amplitude, beat frequency, and chirp 
rate. Once phase-locked, the local Langmuir wave will grow monotonically to large amplitude, then 
saturate with very little ringing. Not all spatial regions will reach precisely the same final amplitude, 
but eventually the wave can grow monotonically more or less everywhere until local saturation ensues, so 
the final variations should be considerably less than in standard approaches. Although plausible, given 
what has been demonstrated in previous analyses of autoresonant phenomena, this expectation should 
also be verified in more realistic simulations. 

We have also neglected all thermal effects, but, in realistic situations, tunneling ionization by and inverse 
bremstrahlung of the drive lasers will typically result in plasma temperatures of T e > O(10 eV). Damping 
of the high phase-velocity plasma waves should remain small, but thermal effects can lead to increased 
particle-trapping and lower thresholds for wave-breaking as well as induce changes in the Langmuir 
dispersion relation, while also resulting in some background of random short-scale density fluctuations 
and electrostatic oscillations which may impact the initial stages of autoresonant phase-locking, which 
assumes a suitably quiescent plasma. This too may be further invesitgated via numerical simulation. 

The suitability of the excited plasma waves for relativistic particle acceleration depends not only on the 
magnitude and uniformity of the peak electric fields fields but even more crucially on the uniformity of the 
phase and phase velocity, and on the degree of phase-locking to the external drive. The Langmuir phase 
velocity v p approximately matches the laser group velocity v g , so first one should consider variability in 
the latter. Diffractive effects can substantially lower the longitudinal group velocity of the laser [19|, but 
once accounted for, additional variation due to the finite bandwidth Slo ~ lu p are typically quite small 

in underdense plasmas, — ~ ^|^, as are variations due to moderate density variations: — ~ Sr— . 
The relative change in the dephasing length (for relativistic electrons with initial velocities very near 
c) is correspondingly small: ~ 

With variations in the group velocity v g expected to be small, phase coherence of the plasma wave will 
depend on how closely v p follows the essentially constant v g — v g of the laser. Particle-in Cell (PIC) 
simulations of Filip et al. 0] suggest that for the RL scheme, the effective phase velocity of the nonlinear 
plasma wave (measured in terms of the progression of the field maximum) can vary appreciably, i.e., 10% 
to 20%, reflecting phase slippage of the Langmuir wave primarily as a result of ponderomotive blowout 
and relativistic detuning, while the plasma wave produced in their non-resonant PBWA scheme exhibits 
substantially less phase slippage. By working within the QSA, in ID geometry without transverse 
density variation, we cannot independently assess any such slippage effects for the present scheme, but 
we anticipate that it will be similarly small by virtue of the autoresoant phase-locking. That is, in the 
resonant RL/TD scheme, plasma inhomogeneities can lead to accelerating buckets that have a changing 
phase, so that electrons do not experience a constant accelerating field. With autoresonance, however, 
the laser can phase lock to a range of densities, creating an accelerating field of uniform phase that is 
everywhere directly related to the local phase of the driving beat wave. 

An appealing feature of both the non-resonant PBWA and autoresonant PBWA is this ostensible ability 
to phase-lock the plasma wave to the beat-wave of the applied drive lasers, with the implied hope that the 
electron injection, whether based on external cathode injection |50l | or internal optical injection |5lLl52| 
can also be phase-locked to the same lasers. Because of its potential importance, this phase-locking 
deserves careful investigation. An obvious worry with the non-resonant PBWA is that entrainment is 
achieved by brute force, in contrast to our self-trapping using adiabatic passage through resonance. 
Because of the intrinsic inefficiency of non-resonant forcing, it must rely on large driving fields and 
denser plasma to achieve large accelerating fields. But increasing u> p and a 2 also increases group- velocity 
dispersion, which can inhibit phase-locking to an externally-known reference phase, and increases the 
growth rates for Raman and self-focusing instabilities which can modulate the laser envelope. One runs 
the risk of turning the non-resonant PBWA into a self-modulated LWFA. That is, such modulation 
can actually enhance the production of plasma waves, but the the plasma response can then become 
entrained not to the initial laser envelope as applied to the plasma, with a prescribed shape and phase, 
but to the envelope after it has undergone some uncontrolled nonlinear evolution and modulation. 
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FIG. 6: Evolving phase difference between the maxima of the laser beat and the longitudinal electric E z (which 
may differ from ^ in the text by a constant). Plot (a) has e = 0.005, a — 0.0001, with the phase slowly evolving 
during autoresonance for < £ < 2700, where E z saturates, (b) uses the CO2 laser parameters of Fig. |3] which 
demonstrates nearly constant phase locking during autoresonant excitation for 40 < £ < 100, permitting good 
control of electron injection. 



One must also carefully distinguish phase-locking from frequency-locking, however much all forms of 
oscillatory entrainment tend to be conflated under the former name. Perfect phase-locking implies 
perfect frequency-locking, and conversely (at least up to some constant but perhaps unknown phase), 
but in the case of only partial or imperfect entrainment, it is possible to achieve good frequency-locking 
without adequate phase-locking, or the converse. Whether the relative error in the phase-matching or 
in frequency-matching between driving and driven oscillation is greater depends on whether the Fourier 
content of fluctuations in the phase is primarily at higher or lower frequencies than the drive frequency 
itself. 

For the purpose of matched particle injection, it is phase-locking which is desired, yet some caution 
is warranted in claims of true phase-locking in either the non-resonant or autoresonant PBWA. In any 
frequency-locked PBWA scheme, the nonlinear frequency of the Langmuir wave may be closely entrained 
to the precisely known drive frequency, but this does not necessarily imply that the absolute phase of the 
Langmuir wave may be precisely known. As a mechanism for phase-locking in the non-resonant PBWA, 
the authors appeal to the claim that an harmonic oscillator, strongly driven off resonance, remains 
synchronous with the driving force, and then suggest that this should extend to nonlinear oscillators. 
But this intuition holds only for damped linear oscillators after the transient is allowed to decay. If a 
linear oscillator, with natural frequency U! n , and negligible damping on the observational time-scales, 
is forced by a constant-amplitude, sinusoidal drive with frequency uJd < w n , then, independent of the 
strength of the drive, the driven oscillator will exhibit persistent, oscillating variations in phase relative 
to the drive phase tpd = ^dt + tpd{0) which are ~ 0(wd/u; n ), i.e., only first-order in their frequency ratio. 
There is no obvious reason to expect in general that nonlincarities will improve matters. 

In the autoresonant case, we may also encounter seemingly good frequency-locking either with compa- 
rably good or disappointingly inferior phase-locking, as shown in Fig. [ft] where we plot, for two of our 
previous examples, the phase lag between the beat-wave of the lasers and the plasma wave, numerically 
estimated by the offset between corresponding relative maxima. Fig.EJa) exhibits oscillating phase vari- 
ation which is slowly varying compared to uj p , but of order < 7r/2 in absolute (not relative) magnitude, 
even within the fully nonlinear regime of autoresonance prior to saturation. On any one shot, an injected 
electron bunch, if sufficiently relativistic (with electron trajectories £ e (i") « £ e (0) ) and sufficiently short 
(bunch length Az e <C cj3 v /c), would experience uniform acceleration and little additional energy spread. 
But the exact phase of the maximal accelerating field could not be reliably predicted from one shot 
to the next, so there would be considerable shot-to-shot variation in average energy gain of the bunch 
and no way to reproducibly exploit the peak field, even if the arrival time of the bunch was perfectly 
correlated with the known beat-phase of the drive lasers in vacuum. 

In contrast, in the second case of Fig.^b), as the oscillator passes through resonance, the plasma wave 
phase tightly locks with the drive phase and remains so until adiabaticity is lost and the excitation 
saturates. In this case, phase-locked injection scheme could reliably place electron bunches near the 
maximum acceleration gradient. It appears that the extent of this phase error is determined by how 
deeply the oscillator phase is trapped in its effective potential, which improves with greater drive strength, 
slower chirp rate, and initial detuning further above resonance. 
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The most obvious drawback to the autoresonant APTR scheme, in comparison to the standard RL/TD 
approach, is the added complication of a chirped drive. But the total chirp required is modest, on the 
order of only co p in order to guarantee, in the face of some uncertainty as to the density, adiabatic passage 
through resonance from above. In typical strongly underdense plasmas this corresponds to only a few 
percent of the laser carrier frequency. In solid-state CPA laser systems, this should be relatively easy 
to achieve, because the laser pulse already undergoes optical stretching/chirping and compression/de- 
chirping in a series of gratings, and the gain bandwidth of the crystals is intrinsically large (~ 20% near 
A ~ 1 /im). One need only adjust the CPA optics so as to only partially re-compress the final chirped, 
amplified pulse, in order to leave some residual chirping. 

In CO2 lasers (or other gas laser systems), the gain bandwidths tend to be narrower, but can effectively 
be increased by operating at sufficiently high pressures to doppler-broaden the rotational lines into an 
overlapping quasi-continuum. If the final intensities and fluences needed are below the damage threshold 
for the required optics, one can simply add a pair of gratings similar to those used in CPA systems in 
order to chirp one of the amplified beams [5j|. If the final intensity is too large, then with somewhat 
more difficulty one might arrange a multi-stage system, where an initial seed pulse is chirped in this 
manner before passing through the final amplifier stage, involving a gas under sufficient pressure to cover 
the bandwidth of the chirped seed. More exotically, one might imagine using nonlinear optical effects in 
a gas or plasma cell to achieve the required frequency shifts. 



VII. CONCLUSION 

We have introduced a straightforward, seemingly minor, modification of the DMG scheme for the chirped- 
pulse PBWA, based on the nonlinear phenomenon of autoresonance with adiabatic passage through res- 
onance (APTR), which nevertheless enjoys certain advantages over previous approaches. Rather than 
starting at the linear resonance and chirping downward at some intermediate rate expected to match, 
on average, the beat frequency to the plasma wave frequency corresponding to the growing value of 
the plasma wave amplitude, we start above resonance and sweep the beat frequency downward past 
the resonance, at any sufficiently slow chirp rate, such that the plasma wave frequency beat frequency 
automatically self-locks to the drive frequency, and the plasma wave amplitude automatically adjusts 
itself consistent with this frequency. This autoresonant excitation achieves higher plasma wave ampli- 
tudes at moderate laser intensities, and, most importantly, appears to be much more robust to inevitable 
uncertainties and variations in plasma and laser parameters. Preliminary analysis has been performed 
within a simplified analytic and numerical model, and wake excitation has been studied using realis- 
tic parameters for TkSapphire and CO2 laser systems. The results are very encouraging, and warrant 
extending investigation to higher-dimensional geometries, more realistic plasma inhomogeneities, and 
self-consistent laser evolution, via numerical solution with fluid and particle-in-cell (PIC) codes. 
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